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Abstract 

Methods are described in the framework of light-cone expansion which al- 
low one to determine the Cabibbo-Kobayashi-Maskawa matrix elements \V u b\ 
and \Vcb\ from measurements of the differential decay rates as a function of 
the scaling variables in the inclusive semileptonic decays of B mesons. By 
these model-independent methods the dominant hadronic uncertainties can 
be avoided and the B — > X u £v decay can be very efficiently differentiated 
from the B — ► X c iv decay, which may lead to precise determinations of \V U \\ 
and | Yd, | . 
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The origins of quark masses, quark flavor mixing and CP violation are among the funda- 
mental issues of particle physics. In the standard model, quark flavor mixing is described by 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix and the phase of it is also responsible 
for the CP violation. The CKM matrix appears to be deeply connected to the origin and 
the values of quark masses, since it results from the diagonalization of the quark mass ma- 
trices [0. The CKM matrix is unpredictable except for unitarity in the standard model, and 
precise determinations of the magnitudes and the phases of its elements and a consistency 
checking in the results are of particular importance. 

The magnitudes of two CKM matrix elements \V u b\ and \V c b\ can be determined from 
inclusive and exclusive semileptonic B meson decays ||. Main theoretical uncertainties in 
these determinations arise from hadronic bound state effects on the underlying weak decays. 
Precise determinations of \V u b\ and \V c b\ require methods which reduce or avoid hadronic 
uncertainties. Especially, improved precision in \V u b\ is an even more pressing need. 

The light-cone expansion provides a solution to the problem of organizing nonpertur- 
bative QCD effects on inclusive semileptonic decays of the B meson in such a way as to 
exploit heaviness of the B meson j§-[7| and allows the selection of observables less affected 
by hadronization and thus better suited for precise determinations of \V u b\ and \ V c b\. In this 
framework the leading nonperturbative effect is described by the distribution function of the 
b quark inside the B meson. Several important properties of the distribution function are 
known in QCD. In particular, the distribution function is normalized exactly to unity as a 
consequence of b quantum number conservation. 

In this paper, we propose to determine \ V u b\ and \ V c b\ by measuring the differential decay 



rates as a function of £ g = [y + yV 2 — q 2 + m 2 )/M in the inclusive semileptonic B meson 

decays B — ► X q £i>, where the invariant variable v is defined by v = q ■ P/M, q stands for the 
momentum transfer to the lepton pair, P and M represent, respectively, the four-momentum 
and the mass of the B meson, and m q is the final quark mass with q = u,c. We will 
show that the differential decay rates dT /d£ q are explicitly proportional to the distribution 
function. Thus the dominant hadronic uncertainties can be avoided with these observables 
either by exploiting the known normalization of the distribution function or through the 
cancellation of the distribution function in the ratio of the differential decay rates. These 
methods are complementary to other methods. The main advantage of our approach is 
that the dominant hadronic uncertainties can be avoided, providing model-independent and 
precise determinations of \V u b\ and \V c b\, and the B — > X u iv decay can be more efficiently 
differentiated from the B — > X c lv decay for determining \ V u b\ than the proposed method to 
measure \V u b\ by the hadronic invariant mass spectrum [BH13I1. 

The decay rate for the inclusive semileptonic B meson decay is given by 

G 2 F \V qb \ 2 fhfK m 

where E is the energy of the decaying B meson, k^ u ) and E^ denote the four-momentum 
and the energy of the charged lepton (antineutrino), respectively. V q b is the element of the 
CKM matrix, which is V c b for the b — > c transition induced decay and V u b for the b — > u 
transition induced decay. Our discussion is made in arbitrary frame of reference throughout 
this work. L^ v is the leptonic tensor: 
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IT = 2{k»k v v + kffi - g^k t ■ k v + ie^^kfk?). 
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Wpv is the hadronic tensor: 



W, 



[IV 



B 



UMJim b), 



(3) 



where j^y) = q{y) r y ll {l — 7s)&(y) is the weak current. The B meson state \B) satisfies the 
standard covariant normalization (B\B) = 2E(27r) 3 S 3 (0). 

The nonperturbative QCD effects on the processes reside in the hadronic tensor. In 
general, the hadronic tensor can be constructed of scalar structure functions by Lorentz 
decomposition, and three structure functions W a (u, q 2 ), a = 1, 2, 3 contribute to the inclusive 
semileptonic B decays under consideration with negligible lepton masses: 
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The inclusive semileptonic decay of the B meson involves large momentum transfer over 
most of the phase space because of heaviness of the decaying B meson. Therefore, the 
integral of Eq. ([3]) is dominated by the light-cone distance in the space-time structure. This 
leads to a factorization of the matrix element in Eq. (|3]) into two parts: one characterizing 
the light-cone space-time singularity and another being the reduced matrix element of the 
bilocal b quark operator 0-0]: 



(B\ 



'jMJi{0j\ \B) = 2{S^ p -ie mu p) [d a A q (y)} (B\b(0)^ (1 - l5 )b(y)\B) 
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where S^ av p = g m g v fi + g^gva — 9iw9a/3 an d is the Pauli- Jordan function for a free 

g-quark of mass m q . The light-cone expansion of the reduced matrix element provides a 
systematic way to organize the nonperturbative effects. The leading term of this expansion 
is given by 



(B\b(0)^(l - 75 )%) \B) y , =0 = 2P? ( X dte** p f{t). 

J o 
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The function /(£) is known as the distribution function of the b quark inside the B meson. 
A similar distribution function has been introduced by the resummation of the heavy quark 
Ml. 



expansion 



The light-cone dominance implies that the structure functions are given by 

W 1 = 2\f(£ + ) + /(£_)] , 
W 2 
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[/(£+) - /(£-)] , 



(7) 
(8) 

(9) 



where 



1S Z 
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Equations include the leading twist contribution; higher-twist contributions are ex- 

pected to be suppressed by powers of q 2 . 

The leading nonperturbative QCD effects on the inclusive semileptonic decays is now 
encoded in the distribution function /(£). The detailed form of the distribution function is 
not known. However, taking y = in Eq.(|6|), b quantum number conservation implies the 
normalization condition 



P»(B 6(0)7^(1-75)6(0) B) = l. 



o 2M 2 ' 

The differential decay rate as a function of the scaling variable £ + for B 
calculated in terms of the distribution function as 

dT W M{l tf55( > , , ) 
d£ + - 48tt3 Z + nt + )*{r q /Z + ) 

+ / dq 2 Ju 2 - q*f(£_ 
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where r q = m q /M, v and £_ are related to £ + and q 2 by 



e + M 2 + q 2 
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The kinematic range of ^ + is r q < t; + < 1. Note that £ + is different kinematic variable for 
B -> X„^z/ and S -> X c £i/, defined in Eq. (|TUp . 

The free quark limit corresponds to /(£) = 5(^ — m&/M). Since the heavy 6 quark 
inside the -B meson is almost on-shell, we expect that the distribution function in reality is 
very close to the delta function, which is supported by the analysis based on the operator 
product expansion and the heavy quark effective theory (HQET) |y||nj, indicating that the 
distribution function peaks sharply around £ = m b /M rs 0.93 0,3- 

The /(£_) term is a consequence of field theory, corresponding to the creation of quark- 
antiquark pairs through the Z-graph. The contribution of the /(£-) term, present as an 
integral in Eq. (|12D , is expected to play less of a role, since (i) in the free quark limit 
the integral vanishes and (ii) the dominant contribution to the /(£-) integral at a given 
£ + resulting from the large £_ region, corresponding to the neighbourhood of the upper 
integration limit for q 2 , is suppressed by v 2 — q 2 . This expectation is confirmed by the 
numerical study discussed below. The spectra dT /d^ + are dominated by the /(£+) term for 
both b — >• c and b — ► u transitions and both spectra have a sharp peak at the same peak 
location as /(£+). In the following we will to good approximation ignore the /(£_) term and 
so the differential decay rates dT/d£ + given by Eq. (|T2|) are proportional to /(£+). 
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In the free quark limit, the tree-level and virtual gluon processes b — > u{c)iv generate 
a spectrum nonvanishing at only one point, i.e., £ + = m&/M, solely on kinematic grounds. 
It is gluon bremsstrahlung and hadronic bound state effects that expand the spectrum over 
the whole range of £ + from m q /M to 1. This unique feature implies that measurements of 
the spectra dT/d^ + would offer the intrinsically most sensitive probe of long-distance strong 
interactions. Indeed, we have shown on the basis of light-cone expansion that the spectra 
dT/d£+ are explicitly proportional to the nonperturbative distribution function. Thus their 
measurements would lead to a direct extraction of the distribution function (see further 
below). On the other hand, the spectra dT /d£ + also provide the most straightforward and 
best way to eliminate the dependence on the distribution function since they are proportional 
to it, so that the dominant hadronic uncertainties can be avoided. Moreover, the B — > X q iv 
spectra dY/d^ + are sharply peaked at £+ = m&/M, since the spectra stemming from the 
tree-level and virtual gluon processes would only concentrate at £ + = m b /M and gluon 
bremsstrahlung and hadronic bound state effects smear the spectra about this point, but 
most of the decay rates remain around £ + = m^/M. This implies that the kinematic cut on 
the b —>■ u scaling variable £ + will make a very efficient discrimination between B — ► X u iv 
and B — > X c lv events, even better than the kinematic cut on the hadronic invariant mass of 
the final state. Consequently, measurements of the decay distribution as a function of the 
b — ► u scaling variable £ + would yield, particularly, a high precision \ V u b\ determination. We 
will discuss these points in detail below. 

For the charmless decay B — > X u £i>, we obtain from Eq. (|T2|) 

IT/ \2ffc \ - 1927l3E 1 gXg -> X ^ Mfrt 



with £ u = {y u + Jv* — ql)/M. We have ignored the up quark mass. To avoid confusion, 
quantities for the b — > u ib — > c) decay are indicated by a subscript u (c); in particular, £ + 
is explicitly denoted as £ u ( c ) for the b —* u (6 — > c) decay. 

By integrating Eq. flI6| ) over £ u and using the normalization condition QTTJ) , one gets rid 
of the distribution function, obtaining 



, 1927r 3 E /-i 1 dT(.B X u &/) 



IKbl " Gp^ i, ^ e! • (17) 

This implies that the known normalization of the distribution function allows a model in- 
dependent determination of \V u b\ from a measurement of the weighted integral of the decay 
spectrum dT(B — > X u £v) / d£ u . 

Likewise, for the charmed decay B — > XJ,v, we obtain from Eq. ( TT2| ) 

it/ Pfft \ - 192w3E 1 dT ( B H X ^ 

6 — g%m 6 Cc^( r c/Cc) ' 1 j 



with £ c = (iSc + ^JVc -Qc+ m-D/M. Integrating Eq. flIB|) over £ c and using the normalization 
condition (|TTD yields 

192^ yi 1 rfr(g -» X c iv) 



G* F M« J rc sc £$(r c /£ c ) d£ c 
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FIG. 1. Contours for £ U)C in the (v,q 2 ) phase space. The dashed (dotted) lines from right to 
left refer to f u(c) = 1,0.9,0.8,0.65,0.5, taking m c = 1.6 GeV. 



which is independent of the distribution function. Thus, a measurement of the weighted 
integral of the charmed differential decay rate as a function of £ c would provide a model 
independent determination of \V c b\. 

The distribution function cancels in the ratio of the differential decay rates at £ u = £ c 
and we obtain from Eqs. ([TH) and (IS) 



ub 



cb 




dT{B -> X u lv)/d£ u 
dT(B -> XJv)/d£ c 



(20) 



This provides a model independent method for determining the ratio | / | . 

Equation fll8| ) shows that a measurement of the decay spectrum dT(B — > X c £v)/d£ c can 
also be used to extract the distribution function directly. We note that the distribution 
function can also be extracted from a measurement of the charmless decay distribution with 
respect to £ u , as shown by Eq. (|16D , though the decay rate is much smaller than the b — > c 
decay rate. 

The quantity £ q (q = u,c) is a combination of two kinematic variables v and q 2 , which 
can be measured through neutrino reconstruction. Except the lower end point £ q = r q 
corresponding to v — and q 2 = 0, different values of v and q 2 correspond to the same 
value of £ g , as shown in Fig. 1, thereby giving the same physical value for the differential 
decay rate dT(B — ► X q £i/)/d£ q . This implies that the spectrum can be well described by 
the light-cone dynamics as long as it can originate from processes with sufficiently large 
momentum transfer. From Fig. 1, we see that this is the case for the dominant spectrum 
over not too small values of £ 9 . Thus departures from the light cone are expected to bring 
little theoretical uncertainty in the methods previously described. Experimentally, one can 
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FIG. 2. The spectra dT(B -> X u £v)/d£, u (left) and dT(B -> X c £v)/d£ c (right) calculated using 
the parametrization of || for the distribution function for a = j3 = 1, a = 0.953, b = 0.00560, and 
m c = 1.6 GeV. The radiative QCD corrections are not included. The absolute scale is arbitrary, 
but the relative scale between the b — > u and b — > c spectra is fixed to |Vub/Vd>| = 0.08. 



take advantage of the multivariate freedom in £ g to choose specific values of v and q 2 in 
phase space so as to discriminate between b —>■ u and b — > c events with the use of kinematic 
cuts and eliminate the charm quark mass uncertainty. If one measures two different pairs 
of (fi,??) and (v 2 , q%) giving the same decay rate dT(B — > X c £u)/d^ Cl this could mean 

V\ + \Jv\ — q\ + Tri 2 = ^2 + ~ <?1 + m ?? which determines m c from experiment, thereby 
eliminating the charm quark mass uncertainty. Note that, as shown in Fig. 1, only the b — > u 
transition is kinematically allowed for ^ c > 1. This cut gives access to a larger portion of 
phase space than the cut directly on v or q 2 . The b — ► u and b — > c events with £ n = ^ c must 
correspond to the different points in phase space, as can be seen in Fig. 1, which renders the 
discrimination between b —* u and b — >• c events easy, in particular when using the method 
of Eq. (g§. 

To verify that the contributions of the /(£-) term in Eq. (|12|) are negligible, we have made 
a numerical study by using the HQET-constrained parametrization of || for the distribution 
function. As expected, the contributions of the /(£_) term to the integrated decay rates in 
Eqs. QP7D and (|T^) are both at the level below 1%, yielding negligible theoretical uncertainties 
in \Vub\ and \V c b\ obtained by these two methods. 

Figure 2 shows the spectra dT(B — > X u £v)/d£ u and dT(B — > X c £v)jd^ c . We see that 
both spectra peak sharply at the same value for £ u and £ c . This feature does not depend 
on the detailed form of the distribution function. However, the shapes of both spectra are 
sensitive to the form of the distribution function. Numerical studies show that including the 
/(£-) term makes almost no difference in the spectra, and the tiny contribution from the 
/(£-) term is meaningful (although still invisible in Fig. 2) only for the spectra very close to 
their two end points, where the spectra become vanishingly small. Therefore, measurements 
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of both £ u and £ c distributions in the neighborhood of maximum, where the /(£-) term is 
negligible, would be appropriate for a reliable determination of \V u b/V c b\ by the method of 
Eq. (pop. The smallness of the /(£-) term has already been found in the parton model M. 

The determination of \V c b\ by the method proposed in this paper may still suffer from 
large theoretical systematic errors. The light-cone expansion parameter is XQ CD /q 2 . For 
B — > X c £i>, the maximum value of q 2 is (M — Mp) 2 with Mp being the D meson mass. It 
means that q 2 is not large enough to neglect the higher order corrections. The AQ CD /g 2 cor- 
rection may be necessary to achieve 5% accuracy for |V^|. Actually the semileptonic b — ► c 
decay rate is dominated by a few exclusive decay modes (D,D* and D^*>ir), which sug- 
gests that the light-cone picture cannot be valid point by point. The theoretical prediction 
obtained in the light-cone expansion is not quite accurate and refers only to the smeared 
spectrum, not point by point. A related problem is the uncertainty in the charm quark 
mass. The charm quark mass is an important source of uncertainty when one measures 
£ c , although, in principle, this uncertainty can be eliminated experimentally by using the 
multivariate freedom in £ c , as discussed previously. Since the spectrum dT(B — > XJ.v)/d£ tC 
has a sharp dependence on £ c , the measurement precision required for the elimination of the 
charm quark mass uncertainty will be challenging experimentally. 

The method proposed in this paper by measurements of the differential decay rate as 
a function of £ u is, on the other hand, theoretically very reliable for a high precision \V u b\ 
determination. The light-cone expansion works much better for the B — > X u iv decay 
because a much larger momentum transfer with the maximum value of q 2 being M 2 can 
occur in the B — > X u tv decay than in the B — > X c tv decay. Many final hadronic states 
contribute to the £ u spectrum above the charm threshold, without any preferential weighting 
towards the low-lying resonance states. Thus, we feel more confident in the method of 
Eq. ([T7|) for determining \V u b\. Applying the kinematic cut v > M — M D , which leads to 
£ u > 1 — Md/M = 0.65 (see Fig. 1), one can discriminate between the inclusive b — > u signal 
and b — > c background. We find that about 99% of the b — > u events pass £ u > 1 — Md/M, 
since the £ M distribution peaks sharply at £ u close to 1, as demonstrated in Fig. 2. This 
discrimination between b — > u and b —>■ c events is even better than the method of the 
hadronic invariant mass spectrum where about 90% of the b —>■ u events survive the cut 
on the hadronic invariant mass Mx < [§-T3|. Thus, b —>■ c background can be very 
efficiently suppressed and improved statistics can be achieved in the measurement of the 
spectrum dT(B — > X u £v)/d£ u , permitting a precise determination of \V u b\. We also anticipate 
that a reliable value for the inclusive charmless semileptonic branching fraction of the B 
meson may be obtained through this measurement. 

The residual hadronic uncertainty in \V u b\ due to higher-order, power-suppressed correc- 
tions is expected to be at the level of 1%. The precision of this determination of \V u b\ will 
mainly depend on its experimental feasibility. The accuracy in v and q 2 , which determine £ u , 
in the experiments seems essential for this method to work. The experimental technique of 
neutrino reconstruction could well make this way of extracting \ V u b\ experimentally feasible. 
If the neutrino can be reconstructed kinematically by inferring its four-momentum from the 
missing energy and missing momentum in each event, then it is possible to measure v and 
q 2 and thus the scaling variable £ u . To reach 10% precision in \V u b\ requires, for instance, a 
measurement of v and y/q 2 with an error around 4%. 

In conclusion, we have proposed new methods for determining \V u b\ and \V c b\ by mea- 



8 



suring the differential decay rates as a function of £ njC in the inclusive semileptonic decays 
B — > X uc iv. These methods are based on the light-cone expansion, which are, in principle, 
model independent as it is in deep inelastic scattering. We have shown that the differential 
decay rates as a function of the scaling variables £ UjC for the inclusive semileptonic B meson 
decays B — > X Ufi lv are proportional to the distribution function, which describes the lead- 
ing nonperturbative contribution, because of the light-cone dominance. This unique feature 
makes these observables ideally suitable for eliminating the dependence on the distribution 
function and thus avoiding the dominant hadronic uncertainties. The integrated decay rates 
are, to a large extent, free of hadronic uncertainties by using the known normalization of 
the distribution function from current conservation, permitting precise determinations of 
\Vub\ and \Vcb\, respectively. Moreover, the distribution function cancels in the ratio of the 
differential decay rates at specific kinematic points, permitting a precise determination of 
The whole methods are then free from model parameters. Since theoretical uncer- 
tainties associated with nonperturbative strong interactions are, to a large extent, avoided 
and perturbative corrections are calculable, these methods could yield very accurate values 
of \V u b\ and \V c b\- There will be smaller residual theoretical uncertainties from higher- twist 
terms and perturbative corrections. 

We wish to emphasize that the method by measurements of the spectrum dT(B — ► 
X u iv) j d^ u is especially promising for a high precision \ V u b\ determination. Theoretically, this 
method is very reliable since the theoretical description based on the light-cone expansion 
works much better for b — > u decays than b — *> c decays. With the dominant hadronic 
uncertainty being avoided, there is a residual hadronic uncertainty in \V u b\ only at the level of 
1%. Experimentally, it gives access to almost all b — > u events, even more than the selection 
against charm background based on the hadronic invariant mass spectrum. Therefore, at 
least potentially, this theoretically sound, clean and experimentally efficient method allows 
for a model-independent determination of \V u b\ with a minimum overall (experimental and 
theoretical) error. 

The decay spectra dT(B — > X u £v) / d£ u and dT(B —>■ X c lv)jd^ c offer, on the other hand, 
the intrinsically most sensitive probe of long-distance strong interactions. Measurements of 
the differential decay rates as a function of £ UjC can be also used to extract the distribution 
function directly. This procedure is crucial for improving the theoretical accuracies on 
the determinations of \V u b\ and \Vd>\ from the charged lepton energy spectra 0|| and the 
hadronic invariant mass spectrum [[0| in inclusive semileptonic B meson decays. A reliable 
value for the inclusive charmless semileptonic branching fraction of the B meson could also 
be obtained by measuring the spectrum dT(B — > X u £v)/d^ u . 
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